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Lambda polarization can be measured through its self-analyzing weak decay, making it an ideal candidate for 
studying spin effects in high energy scatterings. In lepton-nucleon deeply inelastic scatterings (DIS), Lambda 
polarization measurements can probe the polarized parton distribution functions (PDFs) and the polarized frag- 
mentation functions (FFs). One of the most promising facilities for high-energy nuclear physics research is the 
proposed Electron-ion collider in China (EicC). As a next-generation facility, EicC is set to propel our under- 
standings of nuclear physics to new heights. In this article, we study the Lambda production in electron-proton 
collision at EicC energy, in particular Lambda’s reconstruction based on the performance of the designed EicC 
detector. In addition, taking spontaneous transverse polarization as an example, we provide a theoretical pre- 
diction with statistical projection based on one month of EicC data taking, offering valuable insights into future 


research prospects. 
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I. INTRODUCTION 


Spin, as a fundamental property of particles, plays a crucial 
role in the advancement of modern physics. A growing num- 
ber of experimental findings, such as the spontaneous trans- 
verse polarization of A and the proton spin crisis, have made 
it evident that there is much more to be understood about 
spin behaviors in high-energy reactions. The Lambda hy- 
peron (A/A) emerges as an exceptionally powerful tool in 
spin physics, primarily due to its parity-violating weak de- 
cay, which results in a non-uniform angular distribution of its 
products with respect to Lambda’s spin direction [1]. In high 
energy reactions Lambda can be abundantly produced and ef- 
ficiently detected via the decay channel A > pr~ (branching 
ratio is 63.9%). In the A rest frame, the decay protons are 
preferentially emitted along the polarization direction of their 
parent A, with the following angular distribution, 


dN 
T A(1 + aam Pacos), (1) 
where A is the detector acceptance, a, = 0.732+0.014 is the 
weak decay parameter [2], 6* is the angle between proton 
momentum direction and A(A) polarization direction in the 
A rest frame. 

The spontaneous transverse polarization of A was first ob- 
served in 1976 in the unpolarized proton beam scattering on 
Beryllium target [3], by that time the perturbative Quantum 
Chromodynamics (QCD) only predicted a negligible polar- 
ization [4]. This puzzling results triggered a series of theoret- 
ical and phenomenological studies which had been extended 
far beyond A polarization itself. Experimentally, measure- 
ments of A polarization have since been extensively explored 
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in various high-energy processes, encompassing electron- 
positron annihilation [5-7], lepton-nucleon deeply inelastic 
scattering (DIS) [8-10], hadron-hadron scattering [11-13], 
and heavy ion collisions [14-18], yielding invaluable insights 
into numerous aspects of physics. These measurements have 
served diverse purposes, including unraveling the physical 
origins of spontaneous polarization, understanding nucleon 
spin structure, comprehending spin effects in fragmentation 
processes, and exploring extreme conditions of high density 
and high temperature in heavy ion collisions. 


High precision Lambda polarization measurements in the 
proposed electron-ion collider worldwide, provide unique 
Opportunities to study the spin-dependent fragmentation 
functions (FFs) and polarized parton distribution functions 
(PDFs) [19-26]. The Electron-ion collider in China, EicC, 
is the proposed next generation high energy nuclear physics 
facility, which is based on the High Intensity Heavy-ion Ac- 
celerator Facility (HIAF) in Huizhou, China [27, 28]. It is 
conceptually designed to deliver high luminosity electron- 
proton, electron-ion collision, with electron, proton, and 
light ion beams highly polarized. With complementary kine- 
matics coverage to the other electron-ion collider proposals 
worldwide [29-31], the featured physics at EicC includes 3- 
dimensional proton spin structure, nuclear partonic structure, 
exotic hadron states, etc. Lambda polarization measurement 
at EicC is expected to be sensitive to not only the spin de- 
pendent parton distribution functions, but also the spin de- 
pendent fragmentation functions. Potential measurements of 
Lambda transverse polarization and impact studies have been 
performed for US-based EIC which is designed to collide 
electron and proton/ion beams at significantly higher energies 
than EicC [26]. 


In this work, the Lambda production in electron-proton 
scattering under EicC configuration is studied. Based on the 
current conceptual EicC detector design, especially the de- 
sign of tracking subsystem, the reconstruction performance 
for A/A is assessed. In section II, we will describe the simu- 
lation setup including the event generator in use, the detector 


configuration and the corresponding fast simulation proce- 
dure. Performance of A/A reconstruction will be presented 
in section III. In section IV, taking the spontaneous trans- 
verse polarization as an example, the potential statistical pre- 
cision for polarization measurements will be given together 
with theoretical predictions. We will give a brief summary 
and outlook in section V. 


Il. SIMULATION FRAMEWORK 


To simulate the Lambda production in electron-proton scat- 
tering, event generator PYTHIAeRHIC [32], a modified ver- 
sion of PYTHIA6.4.28 [33], is used with the Parton Dis- 
tribution Functions (PDFs) input from LHAPDF [34]. The 
collision energy we choose is the baseline energy outlined 
in the EicC whitepaper [28], 3.5 GeV electron on 20 GeV 
proton. The leading-order diagram for A production in DIS 
process is shown in Fig. 1. The kinematics of the studied 
DIS events are constrained in the following ranges: Bjorken- 
xz 1078 < xp < 1, transferred 4-momentum squared Q? >1 
GeV?, and hadronic invariant mass squared W? > 4 GeV?. 
10 millions of such DIS events are generated for the following 
studies. 


p(P) a x 


Fig. 1. Leading-order diagram for A production in a semi-inclusive 
DIS process. 


At generator level, the average number of A produced per 
DIS event in the above kinematics ranges is about 0.1. In 
the laboratory frame, the momentum and polar angle distribu- 
tions for A and the decay products are shown in Fig. 2. Com- 
paring distributions of daughter proton and pion with of A, it 
can be found that proton carries most of A’s momentum while 
pion only shares a small fraction. A is preferentially produced 
in proton going direction and with a large amount produced 
at very forward angle. Same distributions in Fig. 2 for A are 
similar with slight difference which will be discussed in later 
context. 

In this work, we are mostly interested at the A/A from 
the struck quark fragmentation (current fragmentation re- 
gion). Typically, Feynman-x is an effective variable to sep- 
arate current fragmentation region and target fragmentation 


region. Feynman-z z p is defined as zp = ap) sy, where 
pi) is the A(A) longitudinal momentum in the hadronic 


center of mass frame, and W the hadronic invariant mass. 
Criteria xr > O is expected to suppress the contributions 
from target fragmentation region. The correlation between 
Feynman-x and A/A pseudorapidity 7 are shown in upper 
panels in Fig. 3. Here, pseudorapidity is defined as 7 = 
—In(tan(@/2)), where 0 is the polar angle. Following the 
EicC convention, the positive 7 is along the moving direction 
of proton/ion beam. One can see that A/A with xp < 0 are 
mostly produced at very forward region which will be dis- 
carded in the following simulation and analysis. Considering 
the limited coverage of EicC central detector, |n| < 3 is ap- 
plied for A/A and their daughter particles. Transverse mo- 
mentum, pr, vs. 7 for A and A, with zp > 0 and |n| < 3 are 
shown in the lower panels of Fig. 3. By tracing back the full 
event records in PYTHIA, the origins of such A/ ‘A are shown 
in Fig. 4. At EicC energy, about half of A/A are from decay 
of heavier hyperons. There are also significant contribution 
from beam remnants (di-quarks). In this study, we don’t sep- 
arate different sources for A/A. 


The preliminary conceptual design of EicC detector has 
been described in the white papers [27, 28]. From inner to 
outer, it consists of the vertex/tracking detector, the particle 
identification (PID) system , the calorimeter system, etc. For 
A measurement, the most relevant parts are the tracking and 
PID systems. The latest design of the EicC tracking detec- 
tor is thoroughly described in Ref. [35]. The current design 
of tracking system uses hybrid models. For middle rapid- 
ity (|n| < 1.1), there are 5 layers of silicon and 4 layers of 
Micro-Pattern Gaseous Detectors (MPGD), radially ranging 
from 3.3 cm to 77.5 cm. For |7| > 1.1, the tracking system 
consists of silicon disks followed by large-area Micromegas 
in the forward (proton/nucleus going) direction and all silicon 
disks in the backward (electron-going) direction. For the PID 
system, time of flight detector and Cherenkov detector will be 
used for particle identification at middle and forward rapidity 
respectively. 


For the tracking system, full GEANT4 simulation has been 
performed with the latest design, based on which the resolu- 
tions for primary vertex position, for distance from tracks to 
tracks and from tracks to points, and for track momentum, the 
tracking detector efficiency as function of track pr and 7 are 
given in details in Ref. [35] (Fig.4-9 therein). As well, a fast 
simulation framework is also developed to simulate the de- 
tector responses learned from the GEANT4-base simulation. 
In this work, we follow the same fast simulation procedure 
described in Ref. [35]. The detailed GEANT4 simulation for 
PID system is not available when this work is performed. To 
mimic the particle identification imperfection, a simplified 
“PID smearing” is included in the detector effect fast simu- 
lation. In principal the PID efficiency is correlated with the 
momentum of particles. However, we employ a toy model to 
study the PID effect with a typical PID efficiency of 95% as 
the following. Basically the identified 7, K, or p has 95% 
possibilities to be correct, and 2.5% possibilities to be one of 
other two particles respectively, as described by the following 
matrix: 
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Fig. 2. Momentum (radial) and polar angle (polar) distributions for A and its decay products in the laboratory frame. 


Fig. 3. Upper panels: Feynman-x xp vs. pseudorapidity 7 for A 
(left) and A (right) in the laboratory frame. Only A/A above the red 
line (xr > 0) are kept. Lower panels: transverse momentum pr vs. 
pseudorapdity 7 for A (left) and A (right) with x > 0 and |n| < 3 
(also |7| < 3 for daughter proton and pion). 


T 0.95 0.025 0.025| |7 
K = |0.025 0.95 0.025| |K (2) 
P leaca 0.025 0.025 0.95 P| cones, 


Here 95% of the PID purity is specifically chosen, and a few 
other numbers are also checked for a complete study. 


Il. LAMBDA RECONSTRUCTION 


Similar to the method used in other experiments with track- 
ing detector, A/A reconstruction in this study is based on 
the topological structure of the decay channel with a large 
branching ratio, A —> pr~ and A + pr”. Taking A as ex- 
ample, Fig 5 schematizes the main topological features of its 
production and decay process in a tracking detector. Blue dot 
at the bottom-left represents the ep scattering vertex, named 
as “primary vertex”. A/A is emitted from the primary vertex, 
then moves along the magenta straight dash line and decay 
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Fig. 4. Origins of A and A with xr > 0, |n| < 3 based on PYTHIA 
event record. 


at “VO vertex”. The decay products pr~ (pr™) travel along 
helical lines with opposite bending directions in the magnetic 
field. 

Reconstruction of A/A starts with pairing proton and pion 
tracks with opposite charge. To select the A/A candidates 
and suppress the random backgrounds, the following selec- 
tion variables are considered : 

(1) The distance of closest approach (DCA) of proton and 
pion tracks to the primary vertex. As indicated in Fig. 5, 
DCA, and DCA, from signals should be significantly higher 
than those from background, as parent A/A flies certain dis- 
tances from primary vertex before its decay. 

(2) The distance of closest approach (DCA) between paired 
proton tracks and pion tracks. For A/A signal, this variable 
should be consistent with zero within the track space resolu- 
tion. The decay point (VO vertex) is given by the middle point 
of these two tracks at the closest approach, as indicated by the 
brown triangle in Fig. 5. 

(3) The decay length of A/A candidates, which is the dis- 
tance between the primary vertex and VO vertex. The charac- 
teristic decay length of A hyperon cr is 7.89 cm [2]. 


(4) The angle between A/A candidate momentum jp and 
its trajectory 7 from primary vertex. For the A/A directly 
produced from the primary vertex, its momentum direction 
is supposed to be along its trajectory from primary vertex. 
Correspondingly, cos(7’- p) should be very close to 1. 


Fig. 5. Topology schematic diagram of A production and its decay 
process through A —> pz”. 


To quantitatively determine the selection criteria, the distri- 
butions of proton-pion pairs from pure A/A sample are com- 
pared with the proton-pion pairs from backgrounds, and the 
comparisons are shown in Fig. 6. Based on the comparisons, 
a set of selection criteria are optimized to balance the back- 
ground fraction and the A/A reconstruction efficiency, in or- 
der to keep as many signals as possible while keeping the 
background fraction at a relatively low level. The numerical 
cut conditions are listed in Tab. 1. 


Variables Cut condition 
DCA, > 0.1mm 
DCA, > 0.5mm 

DCA of pr pair < 0.8mm 
Decay length > 1.5mm 
cos (7: p) > 0.95 


Table 1. The summary of topological criteria for A/A reconstruc- 
tion. 


By implementing the aforementioned selection criteria, we 
successfully obtained a clean sample of A/A candidates. The 
invariant mass spectrum of A/A candidates with kinematics 
cuts of zp > 0, |n| < 3, and za > 0.1 (fractional momen- 
tum of A/A is defined as z4 = =e ) is shown in Fig. 7. 
Here the shown histograms are scaled to an integral lumi- 
nosity of 5 fo~+ which is corresponding to about 1 month of 
EicC data taking. With all selection criteria applied, more A 
over A are reconstructed due to the baryon number enhance- 
ment. There is a clean Gaussian signal peak with very limited 
background. The residual background mainly comes from the 
random combinations of oppositely charged particles and par- 
ticle mis-identification. Invariant mass distribution of these 
background is expected to be linear. With typical side-band 
method, the residual background fraction is estimated to be 


about 2.6% for A and 3.0% for A. Signal mass window is set 
to be within 30 width of Gaussian fit, which is (1.106, 1.124) 
GeV/c?. As for the side bands, they are limited to the re- 
gions far from the mass window to avoid fluctuations by sig- 
nals, but not far enough to escape from the signal peak. The 
left side band is (1.083, 1.093) GeV/c? and right side band 
is (1.137, 1.147) GeV/c?. The background under the signal 
peak is estimated as the sum of the two sides-bands normal- 
ized to the signal window. As mentioned in section II, the 
sensitivity of A/A reconstruction to the PID performance is 
assessed by varying the PID "purity" number. For 100% PID 
purity, the residual background fraction is 1.7%, while for a 
90% case, the residual background fraction increases to 3.4%, 
which is still under good control. 

Figure 8 shows the A and A reconstruction efficiency ver- 
sus transverse momentum. The reconstruction efficiency in- 
volves several effects including topological cuts, detector ac- 
ceptance, track efficiency which depends on track pr, track 
n, etc. For A/A with large decay length, the number of track- 
ing detector layers the daughter tracks pass through decreases 
and so does the tracking efficiency. The efficiency at very low 
pr is limited by the detector acceptance due to magnetic field. 
Due to low transverse momenta in the forward region (large 
|n|), the efficiency decreases significantly. As already shown 
in Fig. 3, more A than A are produced at large pseudorapidity 
where the reconstruction efficiency is low, which lead to sig- 
nificantly higher efficiency of A compared A to at pr > 0.5 
GeV/c. When pr is larger than 2 GeV/c, efficiency for A re- 
construction increases and gets closer to A since A at middle 
rapidity start to dominate its production. 


IV. SPONTANEOUS TRANSVERSE POLARIZATION 


In this section, we take A/A spontaneous transverse po- 
larization as an example to explore the physics potentials 
of EicC. The theoretical calculation and statistical projection 
based on our simulation results are described in the following. 

The QCD formalism is used in describing A spontaneous 
transverse polarization P, in the semi- inclusive DIS process, 
e (1) +p(P) > e7 (V) + A(pa, Sa) +X. The Trento con- 
vention [36] is followed in the calculation, where the virtual 
photon moves in the positive z direction and the proton moves 
in the negative z direction, and the differential cross section 
can be written as [37, 38], 


da(S Ara? 2 2 
a( AL) = diad Y 1 a: y {Fou 
dix pdydzad"py yQ? 2(1-e) 22B 
+ |Sarlsin(os, — ba) Fir A e, 


(3) 
where y = 2rgM/Q and Q? = —q?, rp = 2 y= AA 
2 = TI J are Lorentz invariant variables, and SA1, PAL 
are the transverse spin vector and transverse momentum of 
the A hyperon, respectively. F4g = FaB(tB, ZA PAL: Q) 
where the subscripts indicate the polarization of proton and A, 
respectively. yy is the spin-averaged structure function, and 
Fointos, — $a) 

UT 


is the spin-dependent term that contributes to 
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Fig. 6. Distributions of topological variables for A/A signal (red) and background (blue) respectively. 
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Fig. 7. Invariant mass distributions of A and A candidates passing 
all selection criteria. 


the spontaneous transverse polarization. The experimentally 
measured polarization P, is related to the structure functions 
as follows: 


sin(ds, —Pa) 
po (4) 
Fuu 


Within the usual transverse momentum distribution (TMD) 
factorization, at leading twist the structure functions can be 
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Fig. 8. Reconstruction efficiency of A and A as a function of pr 
after all selection criteria applied. 


written as, 


wr 


Fou = [pak (ean, +ki -Pau 


x 5 e+ fia(®B, P}, Q) Di, (za, ki 
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Q), 


rop S = k Pp dki’ (zapi +h — pan) 


BaL 
oe : ZA Ma Pas BL (eB, pè, Q)DEA (za, k2,Q), 


(5) 
where p; and k, denote the transverse momentum of the 
quark relative to the initial proton and the transverse momen- 
tum of A relative to its parent quark, respectively. 
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Fig. 9. The statistical projection with theoretical predictions for A and A polarization in ep collisions at EicC. The size of the projected 


statistical errors are smaller than the marker sizes, thus invisible. 


We parameterize the TMDs using the usual Gaussian form, 
as the product of collinear functions and Gaussian widths: 


eP? / (o3) 


fial(£B, P1: Q) = fig(xB, Q) 


T(pi) 
A 2 A e Fi /(hi) 
Dy, (2a; 13 Q) = Digla Q) Ty (6) 
eo ki / (Mp) 


LA 2. . = LA 
Dirala k1; Q) = Ditg(Za, Q) Ey 


where (p?) = 0.61, (k?) = 0.19 and (M3) = 0.118 
are the corresponding Gaussian widths from Ref. [39, 40]. 
In this analysis, we use the CT18NLO [41] parametriza- 
tion for the collinear PDF while using the DSV [42] and 
AKKO8 [43] parametrizations for the collinear unpolarized 
FF. Both parametrizations describe experimental data but sig- 
nificantly differ, with AKKO8 including substantial isospin 
symmetry violations and DSV parametrization conserving it. 
Additionally, the universality of the polarizing FF DER has 
been proven [44—46]. Similarly, Dit as a modulation of 
DA by an additional collinear function, is also described by 
two different parametrizations, i.e. CLPSW [47] considering 
isospin symmetry and CKT allowing isospin symmetry vi- 
olations [40]. The future EicC experiment provide an ideal 
place to test the A isospin symmetry of FFs. In our study, 
we employ these different parametrizations to calculate the 
polarization observables and perform comparisons. 


Using these parametrizations for the TMDs in Eq.(6), the 
spontaneous transverse polarization of A in Eq.(4) has the an- 


alytic form: 


Pi (TB, ZA, Par, Q) 


— Ly eadig eB, Q)Dihy(za,Q) (kL) + 2% (o) 
Yq Cfig(tB, QDC Q) (M3) + 22 (p2))? 


f2 PAA 
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za Ma 


(7) 
With this expression, we can estimate the magnitude of Pa in 
SIDIS. Taking Q? = 5 GeV’, we plot the P4 as a function 
of pa, in Fig. 9. The results are obtained for different values 
covered by the kinematic range of the future EicC. To get Pa 
dependence on the Feynman variable xr, we parameterize 
zr as a function of Lorentz invariant variables (ap, z4, Q) 
through a kinematic transformation: 


pe —zrQ? | Q24 Qt 


4x? M? 


M[zrB M? + (1— zB)Q?] 


(8) 
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Using Eq.(7,8) and integrating over p, ,, we plot the result 
of the zp and z, dependent P; in Fig. 9. 


The statistical projection of A/A polarization is based on 
an integrated luminosity of 5 fo~!, which is of same size of 
data sample as shown in Fig. 7. The statistical uncertainties 
follows the equation format as ôP = PTE based on the 
polarization extraction procedure. The A and A projected pre- 
cision versus pal, £p, and Zz, are also shown in Fig. 9 to- 
gether with theoretical predictions. The size of the error bars 
are smaller than the marker sizes, thus invisible. Depending 
on the statistics in different bins, the range of the errors is 
from 0.002 to 0.007. 


V. SUMMARY AND OUTLOOK 


EicC is the proposed next generation nuclear physics facil- 
ity which is expected to provide unique opportunities for pre- 
cisely studying the 3-dimensional nucleon structure, the nu- 
clear partonic structure, the exotic hadron states, etc. Lambda 
hyperon serving as a natural final state polarimetry is a pow- 
erful tool for studying nucleon spin structure and spin effect 
in the fragmentation process. Lambda measurements at EicC 
is of special importance and interest. 


Based on a conceptual design of EicC tracking system 
and GEANT4 simulation, we performed a detailed study for 
Lambda production and reconstruction. Also, taking sponta- 
neous transverse polarization as an example, theoretical pre- 
dictions are given as functions of different kinematic vari- 
ables, together with statistical projections with one month’s 
data taking at EicC. We find that measurements with EicC 
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observables will be studied. What could be improved in the 
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